Communication through a Diffusive Medium: 
Coherence and Capacity 
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Coherent wave propagation in disordered media gives rise to many fascinating phenomena as diverse 
as universal conductance fluctuations in mesoscopic metals and speckle patterns in light scattering. 
Here, the theory of electromagnetic wave propagation in diffusive media is combined with information 
theory to show how interference affects the information transmission rate between antenna arrays. 
Nontrivial dependencies of the information capacity on the nature of the antenna arrays are found, 
such as the dimensionality of the arrays and their direction with respect to the local scattering 
medium. This approach provides a physical picture for understanding the importance of scattering 
in the transfer of information through wireless communications. 



The ongoing communications revolution has motivated 
researchers to look for novel ways to transmit information 
(1, 2). One recent development (3, 4) is the suggestion 
that, contrary to long-held beliefs, random scattering of 
microwave or radio signals may enhance the amount of 
information that can be transmitted on a particular chan- 
nel. Prompted by this suggestion, we introduce a realistic 
physical model for a scattering environment and analyt- 
ically evaluate the amount of information that can be 
transmitted between two antenna arrays for a number of 
example cases. On the one hand, this lays a new foun- 
dation for complex microwave signal modeling, an im- 
portant task in a world with ever-increasing demand for 
wireless communication, and, on the other, introduces a 
new arena for physicists to test ideas concerning disor- 
dered media. 

From information theory (5), the capacity of a channel 
between a transmitter and a receiver, that is, the max- 
imum rate of information transfer at a given frequency, 
can be described in terms of the average power of the sig- 
nal S and the noise TV at the receiver: C = log2(l-l-<S'/A^). 
More generally (2), the communication channel connect- 
ing several transmitters and receivers is described by 
a matrix Gia giving the amplitude of the received sig- 
nal a due to transmitter i. The information carried by 
the channel can be characterized using several quanti- 
ties, such as the capacity or mutual information, which 
are typically functionals of the matrix G, which must 
be known in order to predict these quantities. Often G 
cannot be predicted for actual systems, such as wireless 
communication networks or optical fibers, because of the 
complicated scattering and interference of waves that is 
involved. It is crucial, therefore, to develop physical mod- 
els for the signal propagation, as it is only through such 
models that one can understand the real effects of scat- 
tering and interference on the amount of information that 
can be communicated. 

In many cases, only partial information is available for 
prediction; in these situations, one only has a statistical 
description of G. Instead of making assumptions about 
G directly, the usual procedure in information theory, we 
introduce statistical models for the physical environment 



from which we derive the properties of G. The advantage 
of this procedure is that simple physical models can yield 
very non-trivial properties of G. 

Statistical descriptions of the environment have been 
quite successful in the physics of disordered media (6 - 
8). The simplest of these is diffusive propagation. In 
our case of electromagnetic propagation in the context of 
wireless communication, diffusion is known to work well 
in various circumstances (9), and simple extensions seem 
relevant for many others. From a diffusive approach, one 
finds the moments of the distribution of G. These will en- 
able us to calculate information-theoretic quantities (for 
example the capacity) using a replica field theory ap- 
proach to random matrix theory (10). Implicit in this 
approach is the assumption that the full distribution of 
G is sampled, which is realistic in many real world situa- 
tions where the environment is changing. However, when 
the number of antennas is large, many quantities of in- 
terest become strongly peaked around their average and 
this assumption can be relaxed. 

In a statistical description, the scattering of the signal 
is characterized by the mean-free path, £, corresponding 
roughly to the distance between scattering events. When 
(. is large compared to the wavelength A but small com- 
pared to the distance d between the two arrays, the wave 
propagation becomes diffusive (7, 8). This has been an- 
alyzed in the past in the context of electron diffusion in 
metals (6) and light propagation in solids (7, 11). In the 
case of wireless propagation, with signals in the 2 GHz 
region, A ~ 10-15 cm, while (, is of order meters for in- 
doors and tens of meters for outdoors propagation, so 
diffusion is applicable. 

In the diffusive regime A <C ^, to leading order in 
A/£, only the quadratic correlations {G^^G*^) are non- 
negligible and therefore describe the system, where the 
brackets represent an average over realizations of the dis- 
order. Higher cumulants of G are of higher order in X/ L 
Therefore, the distribution of G is Gaussian with zero 



average 



The leading term in (G^^G*^) is eval 



uated by a summation of so-called ladder diagrams (7) 
corresponding to processes in {G^^G*o) where the waves 
from antennas z to a and from j to propagate through 



the scattering medium along identical paths except for 
segments of order (. at each end. 

In several realistic situations discussed below, the cor- 
relations take a particularly simple form: 
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Here Rap and Tij are matrices describing the correlations 
of the receiving and transmitting antennas, respectively, 
and 5" = Tv{GG'^}/nB, is the average power received at 
each of n^ receivers assuming independent signals from 
each of ut transmitters. R and T are normalized such 
that Tr{T} = riy and Tr{i?} = nij. The factorization in 
Eq. |l| of the receiver and transmitter information reflects 
the dominance of the ladder terms: only the segments 
near the antennas contribute a net phase, so only local 
information about the antennas is relevant. Furthermore, 
under general circumstances, T can be expressed in terms 
of the response Xi(i^, e) of antenna i to an incoming plane 
wave in direction k with polarization e and a weight func- 
tion w(k, e) giving the fraction of incident power in that 
direction with that polarization: 
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(2) 



R is defined similarly. 

Having summarized the coherent diffusion results that 
we will need, we now turn to considering information- 
theoretic quantities. While there are many such quan- 
tities that one could study with the diffusive techniques 
outlined above, for concreteness we focus on a problem 
of relevance to wireless communication applications (2 
- 4, 12, 13). We assume that: (a) the receiver knows 
G, but the transmitter does not, (b) the noise at each 
receiving antenna has the same power N and is Gaus- 
sian uncorrelated, and (c) after a long time, the entire 
space of possible G matrices is explored as the environ- 
ment changes, subject to the statistical properties of Eq. 
n^. Under these circumstances, the maximum time aver- 
aged information transfer rate is given by the so-called 
capacity, formally defined in (14), given by (2 - 4) 

C = max (Tr {\og^[In^ + j^ G Q G^]]) (3) 

w 

where /„^ is the tir x n^ unit matrix and Q is the non- 
negative-definite riT X jiT covariancc matrix describing 
the correlations between the signals from the transmitter 
antennas. Here, the angled brackets again indicate an 
average over realizations of G with correlations defined 
by Eq. P. With the normalization of G given above, the 
constraint on the maximum power transmitted yields the 
condition Tr[(5] < ut- For a fixed covariance Q, the max- 
imum information transfer rate is given by the same ex- 
pression only without the maxg in front. Maximization 
over Q then gives the capacity of the channel. 

The assumptions (a)-(c) given above are appropriate 
for certain real world systems, and have been previously 
studied (3, 4)- Assumption (a) is appropriate if the trans- 
mitter occasionaly sends known signals so the receiver 





FIG. 1. Transmitter and receiver antenna arrays in the 
presence of a single scatterer (mirror). Electromagnetic waves 
interfere at both transmitter and receiver antennas. The pres- 
ence of the scatterer introduces an additional path of wave 
propagation between the two arrays and thereby additional 
interference at both ends. As a result, the receiver can resolve 
the individual signals from the two different transmitting an- 
tennas. 



can determine G, and if G changes slowly enough that 
the channel needs to be probed only very occasionally. 
We also assume that the receiver cannot send informa- 
tion about G back to the transmitter. Assumption (c), 
while necessary to make Eq. formally correct, is not too 
important in practice for systems with a large number of 
antennas where the maximum information transfer rate 
for a given Q is strongly peaked about its average value 
(3, 15). 

An important implication of Eq. S, recently pointed 
out in (3, 4)j is that if the environment is sufhcicntly 
multipath, the capacity increases linearly with the num- 
ber of antennas even when the total transmitted power is 
fixed, greatly exceeding the capacity of line-of-sight prop- 
agation. This result may seem surprising since scattering 
may reduce the received signal. However, the scattering 
also produces many independent paths which interfere at 
both the transmitter and receiver. This interference can 
be exploited by placing several antennas at both ends 
in order to increase the capacity. We now illustrate this 
dramatic result with a simple example. 

Consider the case of two antennas at each end and as- 
sume that a single scatterer (a mirror) is at some point 
between the two arrays (Fig. |l|). Here and throughout 
this paper the Fresnel limit will be assumed: the size of 
the transmitter and receiver arrays, is much smaller than 
V^Ad. In the absence of scattering there is only line-of- 
sight propagation. Due to the small effective aperture 
of the receiving array, it is impossible to resolve the dif- 
ferent transmitting antennas. Thus, the capacity is the 
same as if there were only one transmitter and one re- 
ceiver antenna, except that the noise level is reduced by 
a factor of 2 since it is averaged over 2 receivers. Mathe- 
matically, this corresponds to G being a rank one matrix, 
indicating only one channel of communication. Including 
the mirror, the signal received is now the sum of the di- 
rect line-of-sight and the scattered amplitudes. If the 
two incoming waves are at very different angles, then the 



aperture of the receiving array is sufficient to resolve the 
two waves separately. It is thus possible to distinguish 
the two signals originating from the two transmitters. 
Now, G is of rank two so there are two nonzero terms 
(two channels) in Eq. |[ roughly doubling the capacity. 

Foschini and Gans argue (3) that for many scatterers, 
it is possible for G to become of full rank (13) such that 
the capacity is proportional to the number of antennas. 
An increase in capacity can thus be obtained even if S/N 
is reduced by scattering since this latter effect reduces the 
capacity only logarithmically. In order to derive their 
results, they assume that G is a random matrix with 
completely uncorrelated elements. Their situation (fully 
uncorrelated matrix elements) corresponds to R and T 
being identity matrices in Eq. |l| so that the optimal Q 
in Eq. ^ is identity also. This yields, in fact, an upper 
bound on the capacity which can be realized only in very 
special circumstances. 

We now generalize the results of Foschini and Gans 
(3) to situations with nontrivial correlations {R and T 
not identity) which actually occur in real systems. The 
capacity (14) defined by Eq. with such nontrivial cor- 
relations has never been considered previously. These 
correlations (nontrivial T and R) indicate some redun- 
dancy in both the transmitter and receiver arrays thus 
reducing the capacity. 

Surprisingly, we find that the capacity (Eq. ph can 
be evaluated analytically for a large number of anten- 
nas using replica field theory extensively applied in sta- 
tistical physics (10) once T and R are known (that is 
given (Gj^G*^)). The result of this approach is a set 
of algebraic equations for the capacity in terms of the 
eigenvalues of R and T which will be described in de- 
tail elsewhere. In contrast, attempting to evaluate Eq. 
numerically by generating realizations of G with appro- 
priate covariance would be quite difficult since one must 
then maximize over Q. 

We now analyze illustrative situations involving scat- 
tering media. The simplest case is when scalar waves 
are considered and the antennas are deep inside a uni- 
form and isotropic scattering medium with mean free 
path large compared to the size of the arrays. In ad- 
dition, if we assume that the antennas are ideal, making 
a perfect measurement of the field at one point without 
distorting it in any way, then Xj(k) — exp(ifcok • r^) (no 
e appears for scalar waves). Due to isotropy, i/;(k) — 1, 
and thus using Eq. ^, we obtain Tij — /(fco|ri — Vj\) 
and Rap = /(A:o|rQ — r^|) with f{x) — sin(a;)/x where 
fco = 27r/A and r^ are the antenna positions. When the 
antennas are separated by multiples of A/2, all antennas 
of each array on a straight line, there are no correlations, 
and R and T are unit matrices as assumed by Foschini 
and Gans. Any other configuration has nontrivial corre- 
lations and therefore lower capacity. 

An important consequence of the form of the correla- 
tion matrices is the scaling of the capacity with the num- 
ber of antennas. In the limit of a large periodic array of 
antennas, the eigenvectors of T and R are plane- waves 
(by Bloch's theorem) so one can find the eigenvalues an- 
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FIG. 2. Information capacity per antenna (bits per sec- 
ond per Hertz) between identical antenna arrays in a uniform 
isotropic medium as a function of the number of antennas, for 
three representative situations. In all cases the nearest neigh- 
bor distance is A/2, the signal to noise ratio is 100, and we 
have considered scalar waves. Top: Linear arrays of anten- 
nas. In this configuration the antennas are completely uncor- 
related, yielding the maximal capacity per antenna, C /n = 
5.483 bps/Hz/antenna for this signal to noise ratio. Mid- 
dle: Square arrays. The asymptotic value of the capacity per 
antenna is predicted to be 4.43 bps/Hz/antenna. Note the 
slow convergence of the capacity per antenna to its asymp- 
totic value; this effect is due to the slow decay of correlations 
(/(a;)) of antennas at large distances. Bottom: Cubic arrays. 
Here the capacity per antenna decays as n^^''^ up to log cor- 
rections since only the surface layer of antennas contributes 
to the information capacity. 



alytically. When antennas are placed on a one or two 
dimensional array (Fig. 0), the capacity is proportional 
to the antenna number. However, when they are placed 
in a three dimensional array, the capacity scales as n?'^ 
up to a logarithm, not linearly. The model of a scattering 
medium discussed above gives an intuitive interpretation 
of this result: a finite thickness of antennas on the sur- 
face is sufficient to resolve the incoming waves from the 
transmitting antennas; additional antennas in the inte- 
rior are redundant. Therefore the optimal Q is such that 
less power is transmitted from these interior antennas. 

As a second example we consider ideal dipole an- 
tennas and vector electromagnetic waves. The average 
power propagating is the same as for scalar waves, since 
vectorial diffusion (non-zero spin) decays exponentially 
(11). The correlation matrices are now given by using 
Xj(k, e) = (pi • e) exp(i/cok • r^) in Eq. |, where p., is the 
polarization of antenna i. From the resulting expression, 
one finds that if antennas are positioned on a straight line 
with polarization at an angle = 54.73° (cos^ 6* = i) 
with respect to that line, the correlations again vanish 
when antenna separations are multiples of A/2. In this 
case, then, the ideal situation of Foschini and Gans is 
again realized; however, for any other angle, correlations 
play a role even if the antennas are separated by half 
wavelengths. 
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FIG. 3. Maximum information transfer rate per antenna 
(bits per second per Hertz) between two 1-dimensional 50 an- 
tenna arrays with one placed deep inside a half-infinite diffu- 
sive medium (with antenna spacing .5A) and the other above 
it (with spacing 0.6A) at an angle 9 to the surface. Lower 
curve corresponds to the maximum transfer rate for trans- 
mission signal covariance Qij equal to Sij . Upper curve is the 
information capacity, optimizing over Q in Eq. 3. The signal 
to noise ratio is again 100 and scalar waves are treated. Be- 
cause the scattered signal approaches the antennas from only 
one side (the diffusive medium), the correlations in the per- 
pendicular direction cannot become destructive. As a result, 
the correlations can only be minimized when all the antennas 
are positioned at the same height from the surface. Thus, 
the maximum 5.478 bps/Hz/antenna occurs at 6 — 0, when 
the array is parallel to the surface dividing the diffusive re- 
gion from free space. Note that optimization over Q becomes 
more important when 9 increases and there is higher antenna 
redundancy. 



The simplest case in which the geometry of the envi- 
ronment is taken into account is a half-infinite diffusive 
3-dimensional region without any scattering in the other 
half space (Fig. ^). This is a first approximation to mod- 
eling electromagnetic propagation in a city of high-rise 
buildings. When either of the arrays is close to the sur- 
face of the diffusive region, the correlations between an- 
tennas become qualitatively different. Specifically, when 
one antenna array is outside while the other is deep in- 
side the scattering medium, i(;(k) = (k-z)8(k-z) for the 
outside antenna array (for either scalar or vector waves) 
where z is the normal to the interface and Q is the Heav- 
iside step function. Surprisingly, the maximum capacity 
is achieved by pointing the transmitter directly into the 
diffusive medium rather than at the receiver! 

We close by pointing out several simple extensions 
which make these methods applicable to a wide variety 
of realistic situations. First, more realistic modeling of 
the scattering medium could take into account scatterers 
that are not small compared to the wavelength. A con- 
crete example is the scattering off walls inside a building 
where the scattering is highly anisotropic (9). This case 
can be analyzed within our framework; however, the cor- 



relations no longer have the simple factorized form Eq. 

0- . . 

Second, more realistic antennas can be treated by sim- 
ply inserting appropriate functions Xi(k, e) in Eq. g. 
Even for nontrivial antennas, these functions can be ei- 
ther determined numerically or measured experimentally. 
In this way, both practical antenna designs and the ef- 
fects of antenna interactions can be included. 

Finally, we can consider the case where the transmitter 
also has knowledge of G. In this case, the capacity will 
be greater. Here, the averaging over G occurs after op- 
timizing over Q, a process known as "water-filling" (2). 
This case can also be handled within our framework. 
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